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Disordered  Solids  by  Clare  Yu,  Physics  Department,  UC  Irvine,  Irvine,  CA 

Our  research  in  the  last  year  has  concentrated  primarily  on  two  areas:  ( 1 )  the  competition 
between  interactions  and  randomness,  especially  as  seen  in  the  Coulomb  glass:  and  (2) 
strongly  correlated  electrons  with  emphasis  on  the  Hondo  lattire. 

The  competition  between  disorder  and  interactions  is  found  in  a  wide  variety  of  systems. 
For  example,  it  may  be  possible  to  explain  the  low  temperature  properties  of  glasses  with 
a  model  of  randomly  placed  interacting  defects.  As  1  discussed  in  my  previous  report.  I 
performed  a  numerical  simulation  exploring  this  model  which  has  since  been  published  in 
Physical  Review  Letters1  and  the  Phonon  Scattering  Conference  proceedings.*  (Reprints  are 
enclosed.) 

Another  example  is  the  Coulomb  glass  which  is  relevant  to  doped  semiconductors.  In  a 
Coulomb  giass  electrons  interact  with  each  other  via  Coulomb  interactions  and  they  have 
a  choice  of  different  random  places  where  they  can  sit.  At  high  temperatures  the  electrons 
move  around  a  lot,  but  at  low  temperature  they  don’t  change  their  configuration  much  and 
the  system  is  an  insulator.  My  postdoc  (Eric  Grannan)  and  I  have  performed  Monte  Carlo 
simulations  which  clearly  demonstrate  for  the  first  time  that  as  the  temperature  is  lowered, 
there  is  a  spin-glass  like  phase  transition  in  which  the  electrons  freeze  in  place.  We  find 
that  the  transition  temperature  is  about  an  order  of  magnitude  lower  in  energy  than  the 
typical  interaction  energy  scale.  We  attribute  this  to  screening  of  the  Coulomb  interaction. 
At  low  temperature  the  Coulomb  glass  is  an  insulator  with  a  gap  in  its  density  of  states. 
In  our  simulations  we  were  able  to  monitor  the  development  of  this  “Coulomb  gap"  as  the 
temperature  was  lowered.  This  should  be  observable  in  tunneling  experiments.  A  prep-int 
of  this  work  has  been  submitted  to  Physical  Review  Letters.  Our  simulations  were  done  in 
three  dimensions  and  we  are  currently  trying  to  extend  them  to  two  dimensions. 

Nonlinear  effects  in  glasses  are  another  area  of  interest.  At  the  moment,  we  are  trying  to 
understand  some  low  temperature  experiments  being  done  in  Prof.  Doug  Oshoroff's  group 
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at  Stanford  in  which  they  apply  large  electric  fields  to  a  glass  sample  in  a  capacitor  at  low 
temperatures.  If  they  apply  a  large  dc  field  for  a  day  or  so  to  a  glass  sample,  they  ran  watch 
the  dielectric  constant  decrease  slowly  with  time.  Then  if  they  scan  the  dc  field,  they  notice 
a  dip  or  hole  when  they  pass  the  field  where  they  sat  for  a  long  time.  It  seems  from  our 
preliminary  calculations  that  these  features  cannot  be  explain  with  the  standard  model  of 
noninteracting  two  level  systems.  Interactions  appear  to  be  a  necessary  ingredient.  We  can 
qualitatively  reproduce  these  features  in  numerical  Monte  Carlo  simulations  of  an  Is: rig  spin 
glass.  We  are  currently  trying  to  better  understand  the  physics  involved. 

Our  research  has  also  focused  on  strongly  correlated  electrons  which  play  a  key  role  in 
our  understanding  of  a  wide  variety  of  phenomena  su^h  a?  magnetism  and  heavy  ^'rmion 
behavior.  In  order  to  better  understand  the  physics,  experimentalists  often  substitute  low- 
concentrations  of  impurities  into  these  systems.  Dramatic  changes  can  occur  as  a  result. 
For  example,  substituting  in  a  few  percent  of  nonmagnetic  La  ions  into  the  mixed  valence 
compound  CePd3  can  dramatically  increase  the  low  temperature  resistivity.3  But  before  we 
can  understand  the  effect  of  impurities,  we  must  first  understand  the  ordered  lattice.  To  this 
end,  Steve  White  and  I  have  studied  a  one-dimensional  Kondo  lattice.  We  used  a  powerful 
new  renormalization  group  technique  to  calculate  low-lying  energy  levels.  Because  this 
technique  can  calculate  things  that  other  techniques  cannot,  we  have  been  able  to  uncover 
new  physics.  A  preprint  of  this  work  has  been  submitted  to  Physical  Review  Letters.  What 
follows  is  a  brief  description  of  this  wrork. 

The  Kondo  lattice  has  an  f-spin  on  every  site  and  conduction  electrons  that  can  hop 
from  site  to  site  with  a  hopping  matrix  element  t.  If  there  is  one  conduction  electron  on  a 
given  site,  it  interacts  with  the  f-spin  via  spin  exchange  ( J  S;/  -  S,c ).  Spurred  on  by  recent 
experiments  which  have  sparked  interest  in  the  Kondo  insulators,  we  have  examined  the 
half-filled  case  when  the  number  of  conduction  electrons  equals  the  number  of  lattice  sites. 
The  system  is  an  unusual  insulator  in  that  the  gap  to  charge  excitations  is  larger  than  the 
gap  to  spin  excitations.  (In  ordinary  insulators  these  gaps  are  the  same.)  When  the  exchange 


2 


is  much  larger  than  the  hopping  {J  t ),  every  site  has  one  conduction  electron  bound  in 
a  singlet  with  the  local  f-electron.  As  J  decreases,  hopping  becomes  more  important  and 
RKKY  interactions  begin  to  dominate.  RKKY  interactions  promote  antiferromagnetism  and 
cause  f-electron  spins  on  different  sites  to  alternate  in  up-down-up-down  fashion.  We  were 
able  to  study  how  RKKY  interactions  increase  as  J  decreases  by  calculating  the  f-spin-f- 
spin  correlation  function  and  the  staggered  susceptibility.  We  showed  for  the  first  time  that 
the  staggered  susceptibility  diverges  as  J  —0.  We  also  found  a  new  and  unexpected  result 
involving  the  charge  excitations  which  consist  of  a  site  with  two  conduction  electrons  (a 
“particle”)  and  another  site  with  no  conduction  electrons  (a  “hole”).  The  f-electrons  on  the 
particle  and  hob'  sites  f^rm  a  Eaglet.  We  discovered  that  a  particle  and  hole  attract  for  small 
J  ( J  ~  5 1)  but  repel  for  large  J  (J  ~  5<).  Apparently,  for  small  J  the  particle  and  hole  like 
to  sit  next  to  each  other  because  the  RKKY  interaction  has  maximum  attraction  between 


f-electrons  on  nearest  neighbor  sites.  For  large  J,  RKKY  interactions  are  very  weak,  and 
the  particle  and  hole  separate.  The  next  step  in  this  research  will  be  to  see  what  happens 
when  we  remove  the  f-el  con  on  one  site.  This  is  relevant  to  experiments  which  substitute 
nonmagnetic  atoms  in  place  of  magnetic  atoms.  I  also  have  a  student  (Mariana  Guerrero) 


looking  at  the  Anderson  lattice  which  should  be  more  realistic  than  the  Kondo  lattice  in 
describing  heavy  fermion  and  mixed  valence  materials. 

To  summarize,  we  have  studied  (a)  the  competition  between  interactions  and  randomness 
in  the  Coulomb  glass;  and  (b)  strongly  correlated  electrons  in  the  Kondo  insulator.  In  terms 


of  current  and  future  research,  we  are  trying  to  und  ?rstand  nonlinear  effects  in  glasses,  and  . — 

A<3^33<W 

we  plan  to  study  the  effect  of  impurities  on  correlated  electrons.  (j  -  ~ 
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Phase  Transitions  of  Interacting  Clastic  Defects 

Clare  C.  Yu 

Physics  Department.  University  of  California.  Irvine.  California  927/7 
(Received  12  June  1992) 

We  have  invesligated  a  model  of  randomly  placed  defects  vwth  internal  degrees  of  freedom  which  in¬ 
teract  via  elastic  strain  fields.  Monte  Carlo  simulations  and  finite-size  scaling  indicate  that  two  spin- 
glass  phase  transitions  occur:  one  for  the  diagonal  components  of  the  defect  stress  tensor  and  the  other 
for  the  off-diagonal  components.  The  quenched  ground  state  of  the  off-diagonal  components  exhibits  an- 
tiferroelastic  correlations  while  the  diagonal  components  do  not  We  predict  the  fourth-order  elastic  sus¬ 
ceptibilities  associated  with  the  defects  diverge  at  the  transition  temperatures. 

PACS  numbers:  6l.42.+h,  6l.70.Yq.  75.50  Lk.  75.40.Mg 


Defects  in  solids  have  long  been  studied  as  a  model  of 
how  disorder  affects  physical  properties.  In  this  paper  we 
study  the  question  of  whether  a  system  of  interacting  de¬ 
fects  with  internal  degrees  of  freedom  lias  a  phase  transi¬ 
tion  [ 1 1.  In  particular  we  consider  randomly  placed  de¬ 
fects  interacting  with  each  other  via  the  dipolar  elastic 
strain  field  in  three  dimensions.  From  Monte  Carlo  simu¬ 
lations  and  finite-size  scaling,  we  find  that  such  a  system 
undergoes  two  phase  transitions  in  which  degrees  of  free¬ 
dom  freeze  out  as  the  system  is  cooled.  Representing  the 
defects  by  stress  tensors,  we  find  that  the  diagonal  com¬ 
ponents  of  the  stress  tensor  undergo  a  spin-glass  phase 
transition  at  a  slightly  higher  temperature  than  the  off- 
diagonal  components.  “Diagonal"  and  “off-diagona!” 
refer  to  the  projection  oapxaXp  of  the  stress  tensor  on  the 
(cubic)  axes  of  the  underlying  lattice.  In  the  quenched 
ground  state,  the  off-diagonal  components  have  planar 
antiferroelastic  correlations  while  the  diagonal  com¬ 
ponents  have  none.  We  predict  these  two  transitions 
wouid  be  associated  with  diverging  fourth-order  elastic 
susceptibilities.  To  our  knowledge,  this  is  the  first  time 
fourth-order  elastic  susceptibilities  have  been  predicted  to 
diverge. 

We  start  by  considering  defects  that  couple  linearly  to 
the  strain  field. 

H  —  aafi(r) eaP(r) ,  (1) 

where  ea/r)  is  the  symmetric  strain  field  and  aup( r)  is 
the  stress  field  associated  with  the  defects.  The  indices  a 
and  p  range  over  the  real  space  directions  x ,  >\  and  z,  and 
the  sum  over  repeated  indices  is  understood.  As  in  the 
two-level  system  (TLS)  model  of  glasses  [21,  we  assume 
that  the  defects  have  internal  degrees  of  freedom.  Thus 
Cap  can  be  replaced  by  Uap-  S  where  S  is  a  spinlike  TLS 
operator  represented  by  Pauli  matrices.  Vap  is  a  vector  in 
spin  space  and  a  matrix  in  real  space.  The  spin  represen¬ 
tation  is  that  of  the  energy  eigenstates  of  the  two-level 
system.  Sx  and  are  operators  for  transitions  between 
energy  levels,  while  ,Y-  does  not  involve  transitions  and  is 
fsing-like. 

The  defects  interact  with  each  other  via  the  elastic 
strain  field.  For  simplicity  we  just  consider  the  effective 


S:S:  interaction  which  does  not  involve  frequency- 
dependent  effects. 

Using  either  elasticity  theory  [3]  or  second-order  [4,5] 
perturbation  theory  toeliminale  the  strain  field  yields 

—  r')“  —  X  — ^-coslk-  <r  —  r’ ) ] 

\.k  pif 

*  ^p)T}^sciap(r)arS(x')  .  (2) 

where  qip1  mm(iaepk)  +  kee^x)  )/2.  The  sum  over  a.  is  over 
the  longitudinal  and  transverse  phonon  polarizations,  p  is 
the  density  and  r  is  the  speed  of  sound.  ka  is  the  ath 
component  of  the  unit  phonon  wave  vector  and  ep  is  the 
j3th  component  of  the  unit  phonon  polarization  vector 
Summing  over  k  results  in  a  dipolar  interaction  that 
roughly  goes  as  g/r3  where  g—yzlpv 2  and  y  is  the  order 
of  magnitude  of  the  defect  stress.  Taking  y—  1  eV.  we 
estimateg~5x  104  K  A3  for  amorphous  SiO?  (6.7], 

We  simulate  the  system  by  placing  defects  randomly 
on  the  sites  of  a  simple  cubic  lattice.  A  typical  defect 
concentration  is  25%  and  lattice  sizes  were  43,  63,  and  83. 
The  energy  scale  is  given  by  J  ~g/a 3  where  a  is  the  aver¬ 
age  defect-defect  distance.  This  interaction  energy  is  in¬ 
variant  under  block  scaling;  since  the  stress  of  a  block 
scales  as  >//V.  where  /V  is  the  number  of  defects  in  the 
block,  g  scales  as  /V,  and  J  is  invariant  We  mimic  the 
defect  degrees  of  freedom  by  setting  the  magnitude  of  the 
stress  couplings  to  1  eV.  say.  but  allowing  the  sign  of 
each  spatial  stress  component  o„p  to  vary.  In  terms  of 
ro0-  S  this  is  tantamount  to  replacing  S  by  the  magnitude 
of  its  matrix  element  and  treating  each  component  of 
like  a  fixed-length  Ising  spin.  Since  a„p  is  symmetric, 
there  are  six  such  components.  Allowing  the  stress  cou¬ 
plings  to  vary  models  the  ability  of  a  defect  to  respond  to 
fluctuations  in  its  local  strain  field  due  to  temperature 
and  the  fluctuations  of  its  neighbors.  In  the  numerical 
simulations  the  defect  stresses  are  flipped  according  to 
the  standard  heat  bath  Monte  Carlo  algorithm  using  the 
Hamiltonian  (2)  (81. 

The  local  field  of  a  defect  is  determined  bv  summinp 
over  near -neighbor  delects  out  to  a  distance  oi  Rc  ”s/5ao 
where  a0  is  the  lattice  spacing.  We  believe  that  neglect¬ 
ing  further  neighbors  does  not  change  our  results  qualita- 
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tively  for  several  reasons.  First  the  angular  integrals  van¬ 
ish  as  the  stress  due  to  far  away  defects  becomes  more 
isotropic.  Second  if  the  neglected  defects  were  complete¬ 
ly  random  and  one  neglects  the  angular  integrals,  the 
contribution  to  the  local  fields  from  the  remaining  defects 
goes  roughly  as  2gJc(Rcao)~^2  where  c  is  the  fraction 
of  sites  occupied  by  defects.  Since  this  is  of  order  J  for 
r“  25%,  our  calculations  of  the  energy  and  the  specific 
heat  may  be  off  by  a  factor  of  2  or  3.  Finally  simulations 
on  ordered  systems  (c  *■  1 )  indicates  that  summing  over 
32  neighbors  out  to  2a o  vs  56  neighbors  out  to  s/5ao  re¬ 
sults  in  almost  no  change  in  the  ground-state  tnergy  per 
defect  and  the  temperature  at  which  the  specific-heat 
peak  occurs.  The  ground-state  configurations  in  both 
cases  are  antiferroelastic  with  the  off-diagonal  stress  com¬ 
ponents  exhibiting  planar  antiferroelasticity.  An  example 
is  shown  in  Fig.  1.  The  ground  state  is  not  unique,  since 
the  symmetry  of  the  Hamiltonian  allows  up  to  24  degen¬ 
erate  configurations,  though  some  of  these  can  be  identi¬ 
cal. 

Interestingly,  the  tendency  for  planar  antiferroelastici¬ 
ty  survives  in  the  quenched  ground  state  of  the  disordered 
case.  The  metastable  ground  state  is  obtained  by  taking 
the  lowest-energy  configuration  found  at  a  given  tempera¬ 
ture  and  aligning  defect  stresses  along  their  local  fields. 
The  Fourier  transform  of  the  stress  components  quenched 
from  the  lowest  temperatures  is  then  taken  and  averaged 
over  50-100  samples.  For  a  25%  concentration  of  defects 
the  magnitude  of  the  Fourier  transform  of  the  off- 
diagonal  stress  components  \oap(k)\  has  a  peak  at 
kao^Ji  in  the  quenched  ground  state  while  the  diagonal 
components  do  not  show  any  signs  of  ordering.  The  same 
features  are  present  for  10%  concentration  but  the  peak 
height  is  reduced  by  a  factor  of  3  or  4,  indicating  reduced 
correlations  with  increasing  dilution.  Since  the  peak  only 
involves  one  wave  vector  in  both  cases,  this  implies  that 
the  antiferroelastic  correlation  length  is  at  least  as  large 
as  the  system  size  L. 

In  order  to  study  the  phase  transition,  we  define  a  di- 


FIG.  i.  An  example  of  a  ground-state  stress  configuration  in 
the  ordered  case.  Note  the  antiferroelastic  correlations. 


mensionless  order  parameter  ga$  analogous  to  that  used  in 
spin  glasses  by  Bhatt  and  Young  [91, 


Sal 5 


(3) 


where  the  nth  moment  of  the  stress  overlap  qag  is  given  by 


I 


ATfro+D 
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(4) 


Since  the  dipolar  interaction  is  not  random  but  rather  is 
dictated  by  the  relative  positions  of  the  defects,  il.c  aver¬ 
age  I  ]avg  is  over  different  samples  which  have  differ¬ 
ent  placements  of  the  defects.  50-100  samples  are  used 
to  obtain  an  accuracy  of  a  few  percent  in  gag.  i o  is  an  ini¬ 
tial  equilibration  time  and  we  choose  r0”/o  Following 
Bhatt  and  Young  [9],  we  monitored  the  approach  to 
equilibration  by  comparing  the  sample-averaged  overlap 
of  a  set  of  defects  at  different  times  with  the  sample- 
averaged  overlap  of  two  replicas  of  defects  at  the  same 
time.  These  two  overlaps  overestimate  and  underestimate 
the  correlations,  respectively.  Above  the  transition  tem¬ 
perature  these  two  values  converge  to  the  same  value  as 
the  number  of  Monte  Carlo  steps  is  increased,  signaling 
that  equilibrium  has  been  reached.  Below  the  transition 
temperature  convergence  is  hampered  because  the  sym¬ 
metry  of  the  Hamiltonian  allows  degenerate  metastable 
ground  states  with  little  or  no  overlap.  In  this  case  we  ac¬ 
cepted  the  single  replica  value  of  gag  defined  in  (3)  and 
(4)  when  the  fluctuations  that  occurred  in  the  value  of  gae 
were  comparable  to  the  standard  deviation  of  the  sample 
average.  As  a  further  check  of  equilibration,  we  calculat¬ 
ed  the  specific  heal  from  energy  fluctuations  and  checked 
that  values  for  larger  systems  agreed  with  those  for  small¬ 
er  systems  which  were  known  to  have  equilibrated  by  the 
convergence  procedure.  6000  Monte  Carlo  steps  per  de¬ 
fect  were  used  to  achieve  equilibration  for  the  smallest 
sampler  at  the  highest  temperature  while  200000  Monte 
Carlo  steps  were  used  for  the  largest  samples  at  the 
lowest  temperature. 

As  the  size  of  the  system  Ld  — *  the  order  parameter 
gag  varies  between  0  and  1  as  the  temperature  drops 
through  Tr.  According  to  the  finite-size  scaling  ansatz. 
UZ,  is  the  only  relevant  parameter,  where  i  is  the  correla¬ 
tion  length.  This  implies  gap"*gap(L  l/l’<T  —  T,  ))  where 
gag  is  a  scaling  function  and  v  is  the  correlation  length  ex¬ 
ponent.  Figure  2  shows  a  plot  of  gIV  and  gx„  for  l.  “4.  6. 
and  8.  The  curves  cross  at  7,  since  gagiT,  )  is  indepen¬ 
dent  of  L.  Notice  that  Tc  is  different  for  gx,  and  g,v 
since  the  curves  cross  at  different  temperatures.  Fitting 
the  data  by  the  scaling  function  leads  to  7V-/7  *“0.32 
±0.02  and  l/v=  l.0±0.4  for  the  off-diagonal  com¬ 
ponents.  and  T,4-fJ  =0.42  ±  0.02  and  I /v  **0.9  ±0.4  for 
the  diagonal  components  (Fig.  2).  T,  is  lower  for  the 
off -diagonal  components  because  thei'  interactions  are 
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FIG.  2.  (a)  The  x>>  component  of  the  order  parameter  gxy  vs 
the  reduced  temperature  !  ~TU  for  lattice  sizes  L  “4,  6,  and  8. 
The  concentration  c  "•0.25.  The  curves  cross  at  tc-  ”  Tc-/J. 
The  error  bars  at  and  above  72-/7  are  the  larger  of  half  the 
difference  between  the  final  values  of  one-  and  two-replica 
correlations  and  the  standard  deviation  from  sample  averaging. 
The  error  bars  below  Tc-/J  are  the  standard  deviation  from 
single-replica  sample  averaging,  (b)  gxx  vs  t.  The  curves  cross 
at  r<-+**7V+/7.  Notice  that  tc+  >  tc~.  (c)  gxy  fitted  by  the 
finite-size  scaling  formula  with  —0.32  and  l/v  — 1.0.  (d)  gxx 
scaled  with  f<-+  —0.42  and  l/v  — 0.9. 


more  frustrated  than  those  of  the  diagonal  components 
due  to  angular  factors  in  terms  such  as  r,r;rtr/<r,y(r) 
Xffi/(r'),  where  r,  is  the  ith  component  of  the  unit  vector 
connecting  the  two  defects.  To  see  how  the  transitions 
depend  on  input  parameters,  we  doubled  the  longitudinal 
velocity  and  found  that  the  splitting  between  the  transi¬ 
tions  and  l/v  did  not  change  within  the  accuracy 
achieved. 

Each  transition  is  signaled  by  a  divergence  of  the  corre¬ 
sponding  fourth-order  elastic  susceptibility  Xapa»apap  asso¬ 
ciated  with  the  defects.  The  susceptibilities  are  defined  as 
the  coefficients  in  an  expansion  of  the  thermodynamic 
stress  per  defect  in  powers  of  a  small  uniform  external 
strain  field  eap.  This  divergence  is  completely  analogous, 
to  the  divergence  of  the  spin-glass  susceptibility  at  the 
spin-glass  transition  tlOj.  Our  simulations  indicate  that 
the  “spin-glass”  susceptibility  x%fi  ”  A<<?aV  does  indeed 
diverge.  Xa»  is  related  to  fourth-order  elastic  susceptibili¬ 
ties  by  Xapapapap”li)(xap’ ~  f  )  where  p-\/kBT.  We 
have  also  found  the  exponent  p,  which  describes  the 
power-law  decay  of  the  correlation  at  Te,  by  fitting  with 
the  finite-size  scaling  form  x%'  “L  2~nXa0^L  t/v(T 
—  T()).  This  gives  r?-  ■■  -0.2  ±  0.2  for  the  off-diagonal 
components  and  —  +0.4 ± 0.2  for  the  diagonal  com¬ 
ponents  at  c**25%.  When  the  longitudinal  velocity  was 
doubled,  q-  was  closer  to  0.2  but  was  still  0.4.  The 
divergence  of  Xa$  goes  as  l(T ~T()/TC]  'r  where  y 
**(2  - 1?)  v. 

Interacting  defects  have  also  attracted  attention  as  a 
model  for  the  low-temperature  properties  of  glasses.  In 


particular  the  temperature  range  between  3  and  10  K  is  a 
crossover  region  characterized  by  a  plateau  in  the 
thermal  conductivity  and  a  dramatic  drop  by  a  factor  of 
1 02—  1 0 3  in  C/T  with  decreasing  temperature  T  If  we  as¬ 
sume  that  phonons  carry  the  heat  in  this  temperature 
range  ll  1],  then  the  plateau  represents  a  crossover  from  a 
short  mean  free  path  (l—X)  at  high  frequencies  (v >  200 
GHz)  to  a  long  mean  free  path  (/—  I  50a)  at  lower  fre¬ 
quencies  [12].  Attempts  to  explain  this  behavior  have  re¬ 
volved  around  a  drop  with  decreasing  energy  in  the  densi¬ 
ty  of  states  which  strongly  scatter  phonons  [lij.  This 
leads  to  a  drop  in  the  specific  heat  and  a  longer  phonon 
mean  free  path  at  low  energies.  Indeed  neutron  [14]  and 
Raman  [15.16]  scattering  have  seen  evidence  for  this  hole 
at  temperatures  well  ->bove  the  crossover  temperature, 
though  they  have  tended  to  probe  frequencies  and  tem¬ 
peratures  that  are  somewhat  higher  than  those  involved 
in  the  crossover.  A  model  of  interacting  defects  could  ex¬ 
plain  this  since  interactions  tend  to  increase  unperturbed 
energy  splittings  and  produce  a  hole  in  the  density  of 
states.  However,  a  decrease  in  the  density  of  states  can¬ 
not  easily  explain  why  the  ultrasonic  attenuation  at  a 
fixed  frequency  decreases  by  roughly  a  factor  of  30  as  the 
temperature  drops  from  100  to  I  K.  At  low  frequencies 
(viSl  GHz)  this  has  been  attributed  to  relaxation  and 
structural  rearrangement  [2],  At  high  frequencies  [17] 
(v>  100  GHz),  however,  structural  relaxation  is  too  slow 
to  contribute  and  another  explanation  is  needed.  In  par¬ 
ticular  there  is  the  intriguing  possibility  that  the  crossover 
signals  a  phase  transition  in  which  degrees  of  freedom 
freeze  out  as  the  system  is  cooled  [18]. 

We  can  test  this  hypothesis  with  our  model.  Note  that 
for  a  mean  defect-defect  distance  of  a~!5  A  and  y—  1 
cV,  T(~J/3  —  5  K  which  is  the  right  order  of  magnitude 
for  the  crossover  temperature.  We  have  also  calculated 
the  contribution  to  the  specific  heat  from  the  energy  fluc¬ 
tuations  of  the  defects.  As  a  function  of  temperature,  it 
is  a  broad  bump  that  has  a  linear  slope  at  low  tempera¬ 
tures  and  a  maximum  at  T— 7/2.  As  in  spin  glasses  [10], 
Tc±  occurs  at  a  lower  temperature  than  the  maximum. 
The  specific  heat  only  increases  by  a  factor  of  2  from 
7"  "=0.25  J  to  0  5  J  which  is  much  less  than  that  seen  ex¬ 
perimentally  between  3  and  10  K.  Thus  freezing  due  to 
instantaneous  l/r 3  interactions  is  too  mild  to  account  for 
the  low-temperature  crossover  seen  in  glasses.  This  im¬ 
plies  that  the  Hamiltonian  should  include  frequency- 
dependent  interactions  which  involve  defects  making 
transitions  between  energy  levels.  There  are  some  indica¬ 
tions  that  such  terms  would  lead  to  a  greater  drop  in  en¬ 
tropy  upon  freezing  [?]. 

Finally  we  note  that  low-temperature  thermal  expan¬ 
sion  and  ultrasonic  measurements  on  glasses  indicate  that 
defect  stresses  have  a  broad  distribution  centered  very 
close  to  zero  [19).  In  particular  the  ratio  of  the  average 
dilation  stress  per  defect  to  the  rms  value  of  the  stress 
<7>(<t)/3)A<x2)i/2~-  10  ~  \  where  (a2)  is  averaged  over  all 
stress  components.  This  is  naturally  explained  in  our 
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model  by  the  presence  of  disorder.  Averaging  over  1 10 
quenched  ground  states  (/.— 8)  we  find  that  this  ratio 
equals  —0.003  ±0.03  which  is  consistent  with  experi¬ 
ment. 

To  summarize,  we  have  studied  a  system  of  elastically 
interacting  defects.  We  find  that  the  diagonal  and  off- 
diagonal  components  of  the  defect  stress  tensor  undergo 
separate  phase  transitions.  In  the  quenched  ground  state 
the  off-diagonal  components  have  antiferroelastic  correla¬ 
tions.  While  some  of  these  features  may  be  a  result  of  al¬ 
lowing  the  components  of  the  stress  tensor  to  represent 
independent  degrees  of  freedom,  the  divergence  of  the 
fourth-order  elastic  susceptibilities  should  be  valicUfor 
any  elastic  spin-glass  transition  such  as  occurs  in 
KBn  -,KCNt. 
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160  K, 

In  glasses,  the  temperature  range  between  3  and  10  K  is  a  c  rossover  region 
characterized  by  a  plateau  in  the  thermal  conductivity  and  a  dramatic:  drop 
by  a  factor  of  102  —  103  in  C/T  with  decreasing  temperature  T.  where  (’  is  the 
specific  heat  [lj.  In  addition  the  ultrasonic  attenuation  at  a  find  frequency 
decreases  by  roughly  a  factor  of  30  as  the  temperature  drops  from  100  K 
to  1  K.  Below  1  GHz  this  has  been  attributed  to  relaxation  and  structural 
rearrangement  [1].  However,  above  100  GHz,  structural  relaxation  is  too  slow 
to  contribute  and  another  explanation  is  needed  [2j.  In  this  paper  we  explore 
the  intriguing  possibility  that  the  crossover  signals  a  phase  transition  in  which 
degrees  of  freedom  freeze  out  as  the  system  is  cooled  [3]. 

We  model  the  glass  as  a  system  of  randomly  placed  defects  which  interact 
strongly  with  each  other  via  the  elastic  strain  field  [4],  We  start  by  assuming 
that  glasses  contain  defects  that  couple  linearly  to  the  strain  field.  Eliminating 
the  strain  field  yields  [4,5]: 

Heffir  -  r')  =  -  £  — 7  cos(k  •  (r  -  ( 1 ) 

where  oa3{ r)  is  the  stress  field  associated  with  the  defects.  The  indices  a  and 
3  range  over  the  real  space  directions  x,  y  and  z.  =  (k0ejN|  +  kjely1  )/2. 
A  denotes  the  phonon  polarizations  and  e  is  the  polarization  unit  vector,  p  is 
the  density  and  c  is  the  speed  of  sound.  Summing  over  k  results  in  a  dipolar 
interaction  that  goes  as  1/r3  in  addition  to  angular  factors. 

Numerical  simulations  involve  placing  defects  randomly  on  the  sites  of  a 
simple  cubic  lattice.  A  typical  defect  concentration  is  25%  and  lattice  sizes 
were  43,  63  and  $3.  We  mimic  the  internal  degrees  of  freedom  of  the  defect 
by  fixing  the  magnitude  of  the  stress  couplings,  but  allowing  the  sign  of  each 
spatial  stress  component  oa3  to  vary.  We  average  over  50-100  samples;  each 
involving  6000-200.000  Monte  Carlo  steps  per  defect,  depending  on  the  time 
needed  for  equilibration. 

As  we  cool  the  system  down,  it  undergoes  spin  glass-like  phase  transitions 
which  can  be  studied  using  finite  size  scaling  [6],  In  particular  we  find  that 
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t  lit*  diagonal  ami  uir  diagonal  I  «>1 ! !  I  ><  J !  i  t  -  r ! !  s  >,!  !  im  >  S«  •  I  t  '■-I'  ■■  ..'I 

separate  p h as< >  transitions  vv 1 1 It  /  -olf  -  d;av  lUni  ■  /  .-'-e  ■: 

1,1)  2;  11.5  in  hull;  cases.  Em  h  transition  i-  •  •■  •  ■  <•  .:'*•* .  ..•  •  ■  ■' 

1  In’  fm  if  I  li  order  elasi  ic  -uscrpt  liiilit  y  \  ,  .  ,  ,  ,  I  ;,i-  ■  ,<i- '  m  .  •  •  • ..  - 
are  ilflini'il  as  tin*  corlhc  u'lits  in  an  o\pa;is:i  t;  **!  ■  ;:>•  !:.<•:  :iw!i :  :< 

per  lifici  i  in  |)uivi’|s  id  a  '•lliai!  umloini  •  Mei  |  '’la:::  Ins*!  I’  ,  ■  :  ■ 

clear  lnnv  these  ciiver  genres  aifec  t  the  higher  •  '  :*  i  •  • . .  <  - ' :  ■  •••::■•' 
started  with  a  harmonic  Hamiltonian. 

Interestingly.  there  is  a  tendency  lor  planai  ant  i«  ^  ■;.*■ 

olf  diagonal  stress  components  lor  a  25'.  i  concentration  *>! 
nit  tide  of  t  he  funner  transform  *>l  t  in-  uif  diagonal  ri"  •  uiirmi..  ;  :  -  -r  ,  • 

lias  a  peak  at  ka  =  -  while  the  diagonal  rompwiien'  >  *..>  not  >t.u\v  at,-. 
of  ordering  m  the  quern  lied  ground  stale. 

1  he  specific  heal  calculated  from  energy  iim’uatain-  •  A  ‘he  eeam.  >  a 
broad  bump  that  is  linear  at  low  temperature',  with  a  ;i).i\;:;nini  a’  1 
Flie  specific  heat  onlv  decreases  by  a  factor  of  two  horn  !  ..  ...  '■>  1  2 

which  is  much  !«ss  than  that  seen  experimentally  lie’ ween  !  am!  \>l  K  1  lei- 
freezing  due  to  instantaneous  1  interactions  is  too  mild  to  acc  ount  |,,r  tin 
low  temperature  crossover  seen  in  glasses. 

To  summarize,  we  have  studied  a  model  of  glasses  m'.olvin"  lami-uniy 
placed  defects  interacting  via  the  elastic  strain  field.  I'liedefee's  umh'igo  t  v... 
spin  glass  phase  transitions,  1  he  off-diagonal  components  oj  i  In*  .ieiect  :«•" 
tensor  exhibit  antiferroelastic  correlations  in  the  ipienrhed  ground  state  1  ;.<■ 
drop  in  the  specific  heat  is  too  mild  to  account  for  the  crossover  seen  it*  g!as-es 
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